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The influence of the geometry of a thin superconducting 
sample on the penetration of the magnetic field lines and the 
arrangement of vortices are investigated theoretically. We 
compare superconducting disks, squares and triangles with 
the same surface area having nonzero thickness. The cou- 
pled nonlinear Ginzburg-Landau equations are solved self- 
consistently and the important demagnetization efi'ects are 
taken into account. We calculate and compare quantities like 

■ the free energy, the magnetization, the Cooper-pair density, 
\ the magnetic field distribution and the superconducting cur- 

■ rent density for the three geometries. For given vorticity the 
vortex lattice is different for the three geometries, i.e. it tries 

, to adapt to the geometry of the sample. This also influences 
the stability range of the different vortex states. For certain 
magnetic field ranges we found a coexistence of a giant vortex 
placed in the center and single vortices toward the corners of 
the sample. Also the H — T phase diagram is obtained. 
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I. INTRODUCTION 

In mesoscopic samples there is a competition between 
a triangle configuration of the vortex lattice as being the 
lowest energy configuration in bulk material (and films) 
and the boundary which tries to impose its geometry on 
the vortex lattice. For example a circular geometry will 
favour vortices situated on a ring near the boundary and 
only far away from the boundary its influence diminishes 
and the triangular lattice may reappear. Therefore, it 
is expected that different geometries will favour different 
arrangements of vortices and will make certain vortex 
configurations more stable than others. In small systems 
vortices may overlap so strongly that it is more favourable 
to form one big giant vortex. The latter will preferably 
have a circular geometry. As a consequence it is expected 
that the giant to multivortex transition will be strongly 
influenced by the geometry of the boundary as will be 
also the stability of the giant vortex conflguration. 

These issues, the dependence of the stability of the 
giant vortex conflguration and of the different multivor- 
tex configurations on the geometry of the sample will 
be investigated in the present paper. As an example, we 
will compare the most important geometries: the circular 
disk, the square and the triangle. 

Mesoscopic (circular) disks and rings have been the 
most popular in this respect. Experimentally, the mag- 
netization of superconducting disks has been measured 



as a function of the externally applied magnetic field 
. Several transitions between different superconduct- 
ing states were found and the magnetization depends sen- 
sitively on size and temperature. The main part of the 
theoretical studies covered disks and rings P-pX|] of 
zero thickness. In this case one can neglect the magnetic 
field induced by the supercurrents and one assumes that 
the total magnetic field equals the external applied mag- 
netic field, which is uniform. A limited number of studies 
considered disks jrT[-|T4t and rings with finite thick- 
ness. Then, the finite thickness of the disks influences 
the magnetic field profile and it is necessary to take into 
account the demagnetization effects. Often only the (cir- 
cular symmetric) giant vortex states or the superconduct- 
ing/normal transition were investigated. Even in type-I 
superconductors multivortex states in disks |^-P,p^ and 
rings [|l5|,0 were predicted. It was found that if the disk 
or the ring is large enough, the giant vortex nucleates into 
a multivortex state in which the vortices are situated on 
a ring. In a ring geometry, we found that breaking the 
circular symmetry through a non-central location of the 
hole favors the multivortex state [|l5|. This means that 
by changing the geometry, the giant vortex state transits 
into a multivortex state. 

Mesoscopic superconductors with non-circular geome- 
tries have attracted less attention. Moshchalkov et al 
p6| measured the superconducting/normal transition in 
superconducting lines, squares and square rings using re- 
sistance measurements. Bruyndoncx et al [ p^ calculated 
the H — T phase diagram for a square with zero thickness 
in the framework of the linearized Ginzburg-Landau the- 
ory, which is only valid near the superconducting/normal 
boundary. They compared their results with the H — T 
phase boundary obtained from resistance measurements. 
Fomin et al studied square loops with leads attached 
to it and found inhomogeneous Cooper-pair distributions 
in the loop with enhancements near the corners of the 
square loop. Schweigert et al calculated the nu- 

cleation field as a function of the sample area for disks, 
squares and triangles with zero thickness. Jadallah et 
al JTof computed the superconducting/normal transition 
for mesoscopic disks and squares of zero thickness. For 
macroscopic squares, the magnetic field distribution and 
the flux penetration are investigated in detail by exper- 
imental observations using the magneto-optical Faraday 
effect and by first-principles calculations which describe 
the superconductor as a nonlinear anisotropic conduc- 
tor |20|. In the latter case the penetration of the mag- 
netic field occurs continuously. In macroscopic samples 
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the penetration of individual fiuxoids is not so important 
in the theoretical description of the magnetic response 
of the superconductor, but it turns out to be essential 
in our mesoscopic samples. Recently, Aftalion and Du 
pTf studied cylindrical square shaped superconductors 
within the Ginzburg-Landau theory. Chibotaru et al in- 
vestigated the vortex entry and the nuclcation of anti- 
vortices in infinite thin superconducting squares ]2^ ] and 
triangles using the linearized Ginzburg-Landau the- 
ory. Within this linear theory they studied the super- 
conducting/normal transition and they found near this 
transition the nuclcation of multivortices, anti- vortices 
and combination of these two instead of the expected sur- 
face superconductivity. They also calculated the H — T 
phase diagrams for the square and the triangle. Re- 
cently, Bonca and Kabanov |Q studied thin supercon- 
ducting squares using the nonlinear Ginzburg-Landau 
theory in the k ^ oo limit. Within this non-linear the- 
ory they showed that the vortex-anti-vortex configura- 
tion becomes unstable when moving away from the su- 
perconducting/normal transition. 

In the present paper we consider superconductors of 
finite thickness and study also the vortex configurations 
deep inside the superconducting state, i.e. far from the 
superconducting/normal boundary. Our main focus will 
be on the influence of the geometry of the superconductor 
on the vortex configuration and its stability. Our theo- 
retical analysis is based on a full self-consistent numer- 
ical solution of the coupled nonlinear Ginzburg-Landau 
equations for arbitrary value of k. No a priori shape or 
arrangement of the vortex configuration is assumed. The 
magnetic field profile near and in the superconductor is 
obtained self-consistently, and therefore the full demag- 
netization effect is included in our approach. 

The paper is organized as follows: In Sec. II we de- 
scribe the theoretical formalism. Our results are pre- 
sented in Sec. III. We calculate and compare the free en- 
ergy and the magnetization for disks, squares and trian- 
gles with the same surface area. Next, we make a distinc- 
tion between multivortex states and giant vortex states 
and we investigate the influence of the sample geometry 
on the vortex lattice. We also calculate the magnetic 
field range over which the vortex states with vorticity L 
are stable in disks, squares and triangles. The magnetic 
field distribution and the current density are studied and 
the H ~ T phase diagram is obtained. Finally, we sum- 
marized our results in Sec. IV. 



II. THEORETICAL FORMALISM 

In the present paper, we consider thin superconducting 
samples having the same volume but with different geom- 
etry which are immersed in an insulating medium in the 
presence of a perpendicular uniform magnetic field Hq. 
To solve this problem we follow the numerical approach of 
Schweigert and Peeters ||l^]. As for thin disks {d <^ ^, A) 



it is allowed to average the GL equations over the disk 
thickness for samples of arbitrary geometry. Using di- 
mensionless variables and the London gauge divA — 
for the vector potential A, we write the system of GL 
equations in the following form 



AsdA = —S{z) j 
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is the density of superconducting current. The supercon- 
ducting wavcfunction satisfies the boundary conditions 
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and A — ^Hoplt^ far away from the superconductor. 
Here the distance is measured in units of the coherence 
length ^, the vector potential in ch/2e^, and the magnetic 
field in Hc2 = ch/2eS,^ = k^/2Hc- The superconductor is 
placed in the {x,y) plane, the external magnetic field is 
directed along the z axis, and the indices 2D, 3D refer to 
two- and three-dimensional operators, respectively. 

To solve the system of Eqs. l(a,b), we generalized the 
approach of Ref. ]l3[ | for circular disks to superconduc- 
tors with an arbitrary flat geometry. We apply a finite- 
difference representation for the order parameter and the 
vector potential on an uniform Cartesian space grid (x,y), 
with typically 128 x 128 grid points for the area of the 
superconductor, and use the link variable approach [ p5| , 
and an iteration procedure based on the Gauss-Seidel 
technique to find \E'. The vector potential is obtained 
with the fast Fourier transform technique where we set 
A \x\=Rs.\y\=Rs = Hq (x, —y) /2 at the boundary of a box 
with a larger space grid of size typically 4 times the su- 
perconductor area. 

For circular configurations like disks the giant vortex 
state is characterized by the total angular momentum 
L through \I> = ip {p) exp{iL(j)), where p and (j> are the 
cylindrical coordinates. L is the winding number and 
gives the vorticity of the system. Due to the non-linearity 
of the GL equations an arbitrary superconducting state 
is generally a mixture of different angular harmonics L 
even in axially symmetric systems. Nevertheless, we can 
introduce an analog to the total angular momentum L 
which is still a good quantum number and which is in fact 
nothing else then the number of vortices in the system. 

For non-axially symmetric systems there exists no ax- 
ially symmetric giant vortex states and hence the super- 
conducting state is always a mixture of different angular 
harmonics. The vorticity L of a particular superconduct- 
ing sample can be calculated by considering the phase 
(p of the order parameter along a closed loop near the 
boundary of the sample, where the total phase difference 
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is always A</? = Lx 2-k. In non-axially symmetric systems 
three possible vortex states exist: (i) a multivortex state 
which contains separate vortices, (ii) a superconducting 
state which contains one giant vortex in the center, and 
(iii) a state which is a mixture of both: a giant vortex 
in the center which is surrounded by single vortices. The 
giant vortex is not necessary circular symmetric as in the 
case of a circular disk, but it may be deformed due to the 
specific shape of the sample boundary. 

To find the different vortex configurations, which in- 
clude the metastable states, we search for the steady- 
state solutions of Eqs. l(a,b) starting from different 
randomly generated initial conditions. Then we in- 
crease / decrease slowly the magnetic field and recalculate 
each time the exact vortex structure. We do this for each 
vortex configuration in a magnetic field range where the 
number of vortices stays the same. By comparing the 
dimensionless Gibbs free energies of the different vortex 
configurations 



F = V- 



2 M - Ao • J 2d -1*1 
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where integration is performed over the sample volume 
V, and is the vector potential of the uniform mag- 
netic field, we find the ground state. The dimensionless 
magnetization, which is a direct measure of the expelled 
magnetic field from the sample, is defined as 



M = 



47r 



(3) 



where iJo is the applied magnetic field. {H) is the mag- 
netic field averaged over the sample and — rot A . 

The temperature is indirectly included in ^, A, -ffc2j 
whose temperature dependence is given by 
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(4a) 
(4b) 
(4c) 



where Tco is the critical temperature at zero magnetic 
field. We will only explicitly insert temperature if we 
consider the H — T phase diagrams, while the other cal- 
culations are for a certain fixed temperature. Notice fur- 
ther that the Ginzburg-Landau parameter k = A/^ is 
independent of the temperature. 



III. RESULTS 

As a typical example, we consider superconducting 
disks, squares and triangles with the same surface area 
S — 7rl6^^, the same finite thickness d — 0.1^ and the 



same Ginzburg-Landau parameter k = 0.28, which is typ- 
ical for Al [H . Thus the superconducting disk has a radius 
R = 4.0^, the square has a width W = 7.090^ and the 
triangle has a width W = 10.774^. 



A. Free energy and magnetization 

In the first step, we will compare the free energy 
and the magnetization for the three geometries. In 
Figs. |^(a,b) the free energy and the magnetization are 
shown for the disk as a function of the applied magnetic 
field, in Figs. |l](c,d) for the square and in Figs. 0(e,f) for 
the triangle. In Figs. |^(a,c,e) the free energy of the dif- 
ferent giant vortex states is given by solid curves and the 
free energy of the multivortex states by dashed curves. 
The open circles indicate the transition from multivortex 
state to giant vortex state at the transition fields Hmg- 
This transition is of second order Figs. |l|(b,d,f) 

show the magnetization of the different i-states with thin 
curves and the ground state is indicated by the thick 
curve. There exists vortex states with vorticity up to 
L = 11 for the disk and the square and up to 13 for 
the triangle. The superconducting state is destroyed at 
Hc3/Hc2 « 1.95 for the disk, at Hcz/Hc2 ~ 2.0 for the 
square and Hc3/Hc2 ~ 2.5 for the triangle. Thus for 
samples with scharp corners the superconducting/normal 
transition moves to higher field (for fixed surface area) 
[Q . Multivortex states can nucleate in the disk for vor- 
ticity L = 2, 3, 4 and 5 and in the square and the triangle 
for L — 2,3,4,5 and 6. Moreover, for the disk, with in- 
creasing field, the multivortex state always transits to a 
giant vortex state for fixed L, while for the square and 
the triangle some L-states are multivortex states over the 
whole magnetic field range. For the triangle this is the 
case for the states with vorticity L = 3, 4, 5, 6 and for the 
square for the states with L = 4 and 5. This indicates 
that breaking the axial symmetry favors the multivor- 
tex state over the giant vortex state. In some magnetic 
field regions, the vortex states exhibitis a paramagnetic 
response, i.e. —M < 0. This occurs in the disk for 
metastable states with L = 1, 4—8 and in the square 
for metastable states with L = 1, 4, 5. For the triangle 
—M is always positive, i.e. only diamagnetic behaviour 
is observed. 

In Fig. ^(a) we plot the magnetic flux 4>l~>l+i passing 
through the superconductor at the thermodynamic tran- 
sition field Hl^l^i as a function of the vorticity L for the 
three geometries. We scaled this fiux by {L + which 
is the expected flux when no boundary effects are impor- 
tant. The result for the disk is given by circles, for the 
square by squares and for the triangle by triangles. In the 
figure we use open symbols for transitions between two gi- 
ant vortex states, filled symbols for transitions between a 
multivortex and a giant vortex state, and crossed symbols 
for transitions between two multivortices. The symbol for 
the largest L corresponds to the superconducting/normal 
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transition. For the disk we find L ^ L + 1 transitions 
between a giant vortex and a multivortex state for L = 1 
and 2, while the other transitions occur between two gi- 
ant vortex states. For the square we find transitions be- 
tween a multivortex state and a giant vortex state for 
L = 2, 3, 6 and between two multivortex states for L = 4 
and 5. For the triangle, transitions between a giant and 
a multivortex state occur for L = 1, 2 and 6 and be- 
tween two multivortex states for L = 2, 3 and 4. Notice 
further that for L > 6, the value of (j)L^L+i/ {L+l)(j)o be- 
comes almost independent of L for the three geometries: 
4>L^L+il {L + l)(/'o ~ 1.3 for the disk and the square and 
4>L^L+i/ {L + 1)00 ~ 1-4 for the triangle. The magnetic 
flux (j)L^L+i passing through the superconductor at the 
thermodynamic transition field Hl^l+i can be fitted for 
the three geometries by the following formula 

(j)L^L+i _ a + bL 
{L + l)(j)o ^ 1 + cL ' 

with a = 3.10837, b = 0.61792 and c = 0.59825 for the 
disk (solid curve in Fig. |(a)), a = 3.19769, b = 0.60962 
and c = 0.58396 for the square (dashed curve), and a = 
3.58614, b ^ 0.66114 and c = 0.6031 for the triangle 
(dotted curve). 

The amount of flux increase needed to increase the 
vorticity of the ground state from L to i + 1 is plot- 
ted in Fig. 0(b) for the disk (solid curves, open circles), 
the square (dashed curve, open squares) and the trian- 
gle (dotted curve, open triangles). For 3 < L < 11, this 
value is almost independent of L and of the geometry 
and equals A0/(^o — l-lj i-e. it is approximately equal 
to, but still different from, the flux quantum (pQ. Notice 
that > 1 for any L which is a consequence of the 

demagnetization effects. 

B. Multivortex states 

To distinguish whether the superconducting state is a 
multivortex state or a giant vortex state and to determine 
the multivortex to giant vortex state transition field, we 
considered the Cooper-pair density l^*! center center 
of the superconductor ||l^. We can be sure that the su- 
perconducting state is a multivortex state if 1 4*1 center ^ 
for L > 1. The reason is that giant vortices are always in 
the center of the superconductor, and hence 1 4*1 center ~ ^ 
for giant vortex states. Fig. ^ shows the Cooper-pair den- 
sity in the center of the square as a function of the ap- 
plied magnetic field. The Cooper-pair density I'^llenter 
is finite for iJo/i?c2 < Hgm/Hc2 ~ 0.5825 and 0.7825 
for L = 2 and 3, respectively, and l^lcenter — for 
Ho/Hc2 > Hgm/Hc2- For L = A the Cooper-pair density 
in the center differs from zero over the whole magnetic 
field region where the L = 4 state is stable. On the other 
hand l^*] center " ^ ^^^^^ guarantee that the super- 
conducting state is a giant vortex state. For example, 
the multivortex state with L = 5 in a square shows 4 



vortices away from the center situated on the diagonals 
and one vortex in the center, and hence 1^*1 center ~ ^■ 
Therefore, we studied the Cooper-pair density distribu- 
tion in detail. If two vortices are very close to each other, 
then the Cooper-pair density on the axis between these 
two vortices will become very low too, which means that 
the separation between two vortices becomes invisible in 
the contourplots. Therefore we have to define another 
criterion to determine the multivortex to giant vortex 
transition. If the maximum between two minima in the 
Cooper-pair density (i.e. the vortices) is lower than 0.5% 
of the maximum Cooper-pair density 1 5*1 max sam- 
ple, then we will say that the vortices form a giant vortex 
state instead of a multivortex state. With this criterion 
we find that the L = 5 state in a square is always a 
multivortex state and the L — 6 state is a multivortex 
state for Hq/Hc2 < 1-37 and a giant vortex state for 
Ho/H,2 > 1.37. 

How do the multivortex states look like? Figs, ^(a- 
c) show the Cooper-pair density for a multivortex state 
with L = 2,3, and 4 at Hn/Hc2 = 0.42, 0.67, and 0.745, 
respectively. High Cooper-pair density is given by dark 
regions, low Cooper-pair density by light regions. For 
L = 2 the vortices are along the diagonal, for L = 3 the 
vortices are on a triangle, and for L — A they are on 
a square. For L — A only the multivortex state is found 
which is favoured over the giant vortex state. The reason 
is that the square vortex lattice easily fits in the sample. 
Figs. §(d,e) show the phase of the order parameter for 
the multivortex states with L = 5 at Hq/Hc2 — 0.82 
and with L = 6 at Hq/Hc2 = 1.32. Phases near zero 
are given by light regions and phases near 27r by dark 
regions. By going around the superconductor, the phase 
cleary changes 5 times 2tt in Fig. ^(d) and 6 times 27r in 
Fig. ^(e) . For L = 5 there are 4 vortices on a square and 
the fifth vortex is in the center. The latter has clearly 
vorticity one. For L = 6 there are also 4 vortices on a 
square and the other vortices are in the center forming 
one giant vortex with vorticity 2. Thus in this case we 
have the remarkable coexistence of a giant vortex in the 
center with vorticity 2 and 4 clearly separated vortices 
around it. For this case the multivortex to giant vor- 
tex transition field is defined as the field where separate 
vortices appear with decreasing field. This means that 
in Fig. |l| some states are indicated as multivortex, even 
though there exist a (giant) vortex in the center with 
vorticity L > 1. Thus we consider a state no longer as a 
giant vortex state when not all flux of the vortex is con- 
fined in a single connected region. Notice also that not 
only the configuration of the multivortices tries to have 
the same geometry as the sample, but also the giant vor- 
tex geometry depends on the sample geometry. 

We found that for this size of the square sample, the 
states with L > 6 are always in the giant vortex state. 
With increasing L this giant vortex grows and supercon- 
ductivity only occurs in the corners of the square. This 
is illustrated in Figs. |5|(a,b) which shows the Cooper-pair 
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density for the L = 11 state at Hq/Hc2 ~ 1-9 and 1.95. 
It is obvious that further increasing the field pushes the 
superconducting condensate more to the corners. At the 
superconducting/normal transition field Hq/Hc2 ~ 2.0 
the corners becomes normal too. 

Fig. ^ shows the positions of the vortices for the L = 3 
state in a square at applied magnetic fields Hq/Hc2 — 
0.545, 0.62, 0.695, and 0.77. The latter one is just be- 
low the multivortex to giant vortex state transition field 
Hom/Hc2 ~ 0.7825. The solid lines indicate the square 
boundaries. With increasing field the vortices move to- 
wards the center of the square and at Hq m / Hc2 ~ 0.7825 
they combine in the center and form one giant vortex 
with vorticity L = 3. In the multivortex state, one vor- 
tex is always situated on the diagonal of the square, re- 
gardless of the magnetic field. The other two vortices are 
located such that the three vortices form a equilateral tri- 
angle which is centered in the center of the square. Since 
the vortices move to the center with increasing field, the 
width W of the triangular vortex lattice decreases, i.e. 
W = 3.27, 2.89, 2.32, 1.61 at Ho/H^2 = 0.545, 0.62, 

0. 695, 0.77, respectively. 

For the triangle geometry multivortex states nucle- 
ate with vorticity L = 2,3,4,5 and 6 (see Figs. ^e,f)). 
Figs. 0(a-c) show the Cooper-pair density for a multivor- 
tex state with L = 2, 3, and 4 at Ho/Hc2 = 0.495, 0.82, 
and 0.745, respectively. High Cooper-pair density is given 
by dark regions and low Cooper-pair density by light re- 
gions. In the multivortex state with vorticity L = 2 the 
vortices are situated along one of the perpendicular bisec- 
tors of the triangle. In the L = 3 state the vortices are on 
a triangle which easily fits in the sample, while the L = 4 
state consists of 3 vortices on a triangle and the fourth 
vortex is situated in the center. Instead of the square 
configuration as in the case of the square geometry, the 
vortex lattice tries to copy the geometry of the sample, 

1. e. the triangular geometry. For the multivortex states 
with L = 5 and L = 6, the separation of vortices becomes 
invisible in the contourplots of the Cooper-pair density, 
which show one big vortex in the center. The reason is 
that the maximum Cooper-pair density on the axis be- 
tween two vortices is very low. Therefore, we show the 
phase of the order parameter in Figs. 0(d,e) for the mul- 
tivortex states with L = 5 at Hq/Hc2 ~ 1-27 and with 
L = 6 at Hq/Hc2 = 1.345. Phases near zero are given by 
light regions and phases near 27r by dark regions. In both 
cases there is a coexistence of a giant vortex in the cen- 
ter and 3 separated multivortices around it placed in the 
direction of the corners. Taking a loop around the giant 
vortex, the phase changes respectively by 2 and 3 times 
27r, which means that the vorticity of the giant vortex 
is 2 in the case of the L = 5 multivortex state and 3 in 
the case of the L — Q multivortex state. Notice also that 
for L = 6 the geometry of the giant vortex is not axial 
symmetrical, but triangular. States with L > 6 are al- 
ways giant vortex states as we also found for the square. 
With increasing L this giant vortex grows and for large 
vorticities superconductivity only occurs in the corners 



of the triangle. Further increasing the field pushes the 
superconductivity more to the corner untill these cor- 
ners becomes normal too at the superconducting/normal 
transition field. 

For the circular disk we find multivortex states with 
vorticity L = 2, 3, 4 and 5 (see Figs. |l](a,b)). Figs. |(a-d) 
show the Cooper-pair density for the multivortex states 
with vorticity L = 2,3,4,5 at Hq/Hc2 = 0.495, 0.62, 
0.965 and 0.82, respectively. High Cooper-pair density is 
given by dark regions, low by light regions. Multivortex 
states for disks were already studied in previous papers 
(see introduction). Therefore, in this paper we only stress 
that in a disk the multivortices are positioned on a ring, 
which means that also in this case the sample imposes 
its symmetry on the vortex lattice. 

From the study of the Cooper-pair density and the 
phase of the order parameter we learned that: (i) multi- 
vortex states nucleate in disks as well as in squares and 
triangles for several values of the vorticity L, and (ii) the 
vortex lattices tries to have the same geometry as the 
sample. 



C. Magnetic field distribution - Demagnetization 
effects 

Since we studied samples with finite thicknesses, de- 
magnetization effects are important and therefore we had 
to solve for the magnetic field distribution around the 
sample. In this section we will describe the magnetic 
field distribution for the square. The results for the disk 
and the triangle are analogous. 

In Figs, ^(a-f) the magnetic field distribution is shown 
for the square geometry for the state with vorticity L — 2 
at Ho/Hc2 = 0.42, 0.52 and 0.62 (see open circles in 
Fig. |), and with vorticity L = 3 at Ho/Ha2 = 0.62, 0.72 
and 0.82, respectively. High magnetic field is given by 
dark regions and low by light regions. The magnetic field 
is clearly non uniform in and around the sample. The 
dark spots in the square are the vortices and the dark re- 
gions near the sample surface are due to the compression 
of the magnetic lines when they are forced to go around 
the sample. These regions are responsible for the demag- 
netization effects. It is clear that with increasing external 
field and fixed number of vortices, the demagnetization 
effects are more pronounced, because the superconductor 
has to expel more magnetic field. In Figs. ||(a,b,d,e) the 
superconducting state is a multivortex state and the sep- 
arated vortices are clearly visible, while in Figs. |^(c,f) 
where Ho/Hc2 > Hmg/Hc2 ~ 0.5825 and 0.7825 for 
L — 2 and 3, respectively, there is one giant vortex in 
the center. With increasing field the vortices move to- 
wards the center and at Hq = Hmg they combine to one 
giant vortex state. Notice that the giant vortex state is 
not necessarily axial symmetric as was in the case of the 
disk. 

In Figs, l^(a-d) the magnetic field distribution is 
shown for the square geometry for the state with vor- 
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ticity L = 4 at the magnetic fields indicated by the open 
circles in Fig. | i.e. Hq/H^2 = 0.72, 0.82, 0.92 and 1.02, 
respectively. Now, there is no transition from multivortex 
to giant vortex state and the four vortices are cleary visi- 
ble as the dark spots. Notice that from the magnetic field 
distribution one clearly observes that the vortex lattice is 
a square lattice, i.e. the lattice geometry is the same as 
the sample geometry, and that the vortices move towards 
the center with increasing field. 

For the multivortex states with higher vorticity, the 
separated multivortices are not visible anymore in the 
contour plots of the magnetic field distribution. The 
problem is the same as for the contour plots of the 
Cooper-pair density, i.e. the vortices are too close to 
each other and the spots corresponding to high magnetic 
fields are overlapping in the picture. For high vorticity 
and high external fields, the total magnetic field appre- 
ciably differs from the externally applied field only in 
the corners of the square. Figs. |l|(a,b) shows the mag- 
netic field distribution for the same configuration as in 
Figs. ||(a,b), i.e. the L = 11 state at Hq/Hc2 ~ 1-9 and 
1.95, respectively. A local decrease in magnetic field is 
given by the light regions, an increase by the dark regions. 
In both pictures the magnetic field is only substantially 
expelled in the corners and consequently only near the 
corners there is a higher density of magnetic field lines 
at the outside of the square. Further increasing the field 
destroys the superconductivity and, thus, the total field 
becomes equal to the external one over the whole sample. 

Next, we investigate the dependence of the magnetic 
field on z. Figs, ^(a-f) show the magnetic field distribu- 
tion for the L — A state in a square for different values of 
z, z/^ = 0.0, 0.1, 0.3, 0.6, 1.0 and 10.0, respectively. The 
applied magnetic field is Hq/Hc2 = 0.77. High magnetic 
field is given by dark regions, low magnetic field by light 
regions. In the plane of the superconductor, i.e. z — 0, 
the magnetic field which penetrates the superconductor 
is either compressed into multivortices or expelled to the 
outside of the sample. Therefore, the four dark spots in 
Fig. ^(a) indicate the vortices and the light regions to- 
wards the sample boundary are due to the expulsion of 
the magnetic field towards the outside of the supercon- 
ductor. As a consequence, the magnetic field increases in 
a small strip near the sample boundary. With increasing 
z and 1 2; I > d/2, the magnetic field will still be influenced 
by the superconductor. The demagnetization effects de- 
crease with increasing z and the compression of the mag- 
netic field lines into vortices becomes smaller. Therefore, 
the vortices and the expulsion of the field will become 
less pronounced with increasing z. At z = O.li^, the vor- 
tices and the results of the magnetic field expulsion are 
still visible by the dark and light regions (see Fig. p^(b)). 
In Fig. |l^(c) and (d), at z — 0.3^ and 0.6^, respectively, 
the contrast in the picture decreases which means that 
the influence of the superconductor, i.e. the compression 
and expulsion of the magnetic field lines, decreases. At 
z = 1.0^ the magnetic field just shghtly decreases right 
above the superconductor compared to the external field 



(see Fig. |T^(e)). At z = 10.0^ the total magnetic field 
is homogeneous. Far away from the superconductor, the 
magnetic field is not influenced by the superconductor 
and equals the external fleld. This is clearly shown in 
Fig. |l|(f). 

D. Superconducting current density 

When a superconducting sample is placed in an exter- 
nal magnetic field, the magnetic field is expelled from the 
superconductor due to screening currents near the sam- 
ple boundary. The direction of the screening currents is 
such that the corresponding magnetic field is opposite to 
the external one, which leads to a lower total field in the 
superconductor. Magnetic field penetrating the super- 
conductor creates currents flowing in the other direction 
than the screening currents. The competition between 
these currents and the screening currents results in the 
existence of vortices. 

Figs. |l^(a-d) show vector plots of the supercurrent 
in the superconducting square for the L = I state at 
Ho/Hc2 = 0.27, the L = 2 state at Hq/Hc2 = 0.42, the 
L — 3 state at Hq/Hc2 = 0.67 and the L = 4 state at 
Ho/Hc2 = 0.745, respectively. Figs. [l3|(e-h) show the 
corresponding contour plots of the phase of the order 
parameter. Phases near 2tt are given by dark regions, 
phases near zero by light regions. From the phase of the 
order parameter one can easily determine the number 
and the positions of the vortices. In Fig. |l^(a) it is clear 
that the screening currents near the sample boundary 
flow clockwise and the currents around the vortex in the 
center counterclockwise. In Figs. |l3|(b,c,d) there are cur- 
rents flowing counterclockwise around 2, 3 and 4 vortices, 
respectively. Around one vortex, the size of the current, 
indicated by the length of the arrows in Figs. |l3|(a-d), is 
not the same for every angle. In Fig. |l^(b) it is clear 
that in the region between the two vortices the currents 
around these two vortices cancel each other out. Also, 
in the case of L = 3 and L = 4 the currents around the 
different vortices cancel out each other in the center of 
the sample (see Fig. 0|(c,d)). 

From Figs. |l^(a-d)one expects anti-vortices towards 
the corners, because there are some spots where the cur- 
rents flow in clockwise direction. That these are not re- 
ally anti-vortices can be seen from the phase of the or- 
der parameter (Figs. |l^(e-h)). By going around an anti- 
vortex, the phase changes with —2tt and this is clearly 
not the case here. Moreover, the Cooper-pair density, 
shown in Figs, ^(a-c), is not zero at these positions. 

Next, we investigate the superconducting current den- 
sity in the triangular sample. Figs. |lj(a-c) show vector 
plots of the supercurrent in the superconducting triangle 
for the L = 1 state at Hq/Hc2 = 0.27, the L = 2 state at 
Ho/Hc2 = 0.495, and the L = 3 state at Ho/Hc2 = 0.82, 
respectively. Figs. |lj(d-f) show the contour plots of 
the corresponding phase of the order parameter. Phases 
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near 2tt are given by dark regions, phases near zero by 
light regions. The behavior of the supercurrent in tri- 
angular samples is similar to the one in square samples. 
The screening currents flow clockwise and the current 
around the vortices in the opposite direction. The cur- 
rents around different vortices cancel each other in the re- 
gion between them. Towards the corners, there are some 
spots where the current flows in a clockwise direction, 
but these spots are not anti-vortices. This can be seen 
from the phase of the order parameter (see Figs. |l^(d-f)) 
and from the Cooper-pair density (see Figs. |^(a-b)). 

E. Stability of the vortex states 

Not only the stability region of the multivortex states 
with respect to the giant vortex states depends on the 
sample geometry, but also the stability of each individ- 
ual superconducting state is sensitive to the geometry. 
In Fig. |l^ we show the magnetic field range AH over 
which the vortex state with vorticity L is stable, i.e. 

AH = Hexpulsion " HpenetraUon (sCC alsO FigS. |l6|(a,b)), 

as a function of the vorticity L, for L < 6 and in the 
inset for L > 6. For the disk the result is shown by the 
open circles, for the square by the open squares and for 
the triangle by the open triangles where the curves are 
guides to the eye. For the circular disk the stability re- 
gion AH/Hc2 uniformly decreases with increasing L with 
a slight dip at L — 2,3. The square and the triangle ex- 
hibit a peak structure in the region L < 5. For the square 
we find that the state with L = 4 is stable over a larger 
magnetic field region than the state with vorticity L = 3, 
which is a consequence of the fact that the vortex lattice 
tries to keep the same geometry as the sample. For the 
triangle we find a peak at L = 3 and a dip at L = 2 
for the same reason. Notice that: i) the peak structure 
is more pronounced for structures which fit more closely 
the triangular Abrikosov lattice; ii) for L > 4 no clear 
peaks are found; iii) the vortex states in the square and 
circle geometry have almost the same stability range for 
L < 2 and L > 6; iv) for L > 4 the stability range for 
the vortex state in the triangular geometry becomes sub- 
stantially smaller than for the other two geometries which 
have less sharp corners. Thus sharp corners decrease the 
stability range of the vortex states, i.e. it is easier for 
flux to enter the system; and v) Fig. ^, which shows the 
extra flux needed to increase the vorticity by one unit, 
contains complementary information to Fig. 

F. H-T phase diagram 

Untill now, all our calculations were done for fixed 
temperature T. Now we will include temperature and 
our lateral dimensions and fields will be expressed in the 
zero temperature results ^(0) and Hc2{0), respectively. 



Temperature will be expressed in units of the zero mag- 
netic field criticial temperature Tcq. We take the sur- 
face area of our samples S — 167r^^(0) and the thickness 
0.1^(0). 

The H — T phase diagram is shown in the inset of 
Fig. ^ for the disk (solid curves), the square (dashed 
curves) and the triangle (dash-dotted curves) for the 
states with vorticity L = and L = f, thus for low 
fields and temperatures close to Tco- The thick curves are 
the superconducting/normal transitions and the thinner 
curves indicate the expulsion and the penetration fields, 
i.e. the boundaries of the stabilitiy region of the state 
with vorticity L. The lower thin curves show the transi- 
tion from the state with vorticity L = 1 to L = with de- 
creasing field (expulsion) and the upper thin curves show 
the transition from L — Q to L — 1 with increasing field 
(penetration) . Fig. ^ shows the H — T phase diagram 
for higher fields. For the sake of clarity only the super- 
conducting/normal transition iJc3 is shown as a function 
of temperature. The black dots indicate the transitions 
fields between the different L-states. For every (fixed) 
temperature the superconducting/normal transition field 
is highest for the triangle and lowest for the disk. For 
every (fixed) magnetic field, the critical temperature is 
highest for the triangle and lowest for the disk. This 
means that with sharper corners, the critical temperature 
and critical field are enhanced due to an enhanced surface 
superconductivity. These results are in good agreement 
with the phase diagrams found in Ref. [OJ23J23]. 

IV. CONCLUSIONS 

We investigated theoretically the influences of the ge- 
ometry of thin superconducting samples on the vortex 
configuration. Therefore, we considered superconducting 
disks, squares and triangles with the same surface area 
S = 7rf6^^ and the same thickness d = 0.1^ for k, = 0.28. 
For these three geometries we calculated the free energy 
and the magnetization of the different giant and multi- 
vortex states as a function of the applied magnetic field 
and we indicated the multivortex to giant vortex transi- 
tions for fixed vorticity L. Multivortex states were found 
for disks as well as for squares and triangles for several 
values of the vorticity. For given i, the vortex lattice 
was different in the three geometries due to the fact that 
it tries to adapt to the geometry of the sample. This in- 
fluences considerably the stability range of the different 
vortex states. For squares and triangles we found mag- 
netic field regions where there is a coexistence between 
a giant vortex state in the center and several separated 
vortices in the direction of the sample corners. Near the 
superconducting/normal transition we do not find multi- 
vortices, anti-vortices or a combination of them, but we 
flnd surface superconductivity. We studied the magnetic 
held distribution across the superconductor and around 
the superconductor which clearly shows the demagneti- 
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zation effects, which are very important for samples of 
finite thickness. The vector pfots of the superconducting 
current showed spots where the current flows in clock- 
wise direction. From the phase of the order parameter 
and the Cooper-pair density we conclude that these spots 
are not anti- vortices. We also investigated the stability 
of the vortex states with vorticity L by calculating the 
magnetic field range over which the vortex states with 
vorticity L are stable. We found that this stability range 
sensitively depends on the sample geometry. As a func- 
tion of L we found enhanced stability for the triangle for 
L = 3 and for the square for L — A. In the last section, we 
also included temperature by calculating d. H — T phase 
diagram for the disk, the square and the triangle. With 
sharper sample corners, we found that for fixed temper- 
ature, the superconducting/normal transition field H^^ 
moves to higher fields, and for fixed field, the critical 
temperature increases. 
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FIG. 1. The free energy and the magnetization as a func- 
tion of the applied magnetic field for the disk, the square, 
and the triangle with the same surface area S = 7rl6^^ and 
thickness d = 0.1^ for k = 0.28. (a,c,e) The free energy 
of the giant vortex state (solid curves) and the multivortex 
states (dashed curves) and the multivortex to giant vortex 
transition fields (open circles). (b,d,f) The magnetization of 
the L-states (solid curves) and the ground state (thick solid 
curve) . 
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FIG. 2. (a) The flux corresponding to the thermo- 

dynamic transition field Hl^l+i as a function of the vorticity 

L for the disk (circles), for the square (squares) and for the 
triangle (triangles). The curves are the fitted results, (b) The 
fiux A(f> corresponding to the magnetic field range needed to 
increase the vorticity of the ground state from L to L + 1 as 
a function of the vorticity L for the disk (solid curves, open 
circles), the square (dashed curve, open squares) and the tri- 
angle (dotted curve, open triangles). 

FIG. 3. The Cooper-pair density in the center of the square 
as a function of the applied magnetic field for the states with 
vorticity L = 2, 3 and 4. 

FIG. 4. (a-c) The Cooper-pair density for a multivortex 
state in a square with L = 2, 3, and 4 at Ho/Hc2 = 0.42, 0.67, 
and 0.745, respectively. High Cooper-pair density is given 
by dark regions, low Cooper-pair density by light regions. 
(d,e) The phase of the order parameter for the multivortex 
states with L = 5 at Ho/Hc2 = 0.82 and with L = 6 at 
Ho/Hc2 = 1.32. Phases near zero are given by light regions, 
phases near 27r by dark regions. 

FIG. 5. The Cooper-pair density for the L = 11 state in a 
square at Ho/Hc2 « 1.9 (a) and 1.95 (b). 

FIG. 6. The multivortex positions of the L = 3 state in 
a square at applied magnetic fields Ho/Hc2 = 0.545, 0.62, 
0.695, and 0.77. 

FIG. 7. (a-c) The Cooper-pair density for the nmltivortcx 
states in a triangle with L = 2, 3, and 4 at Ho/Hc2 = 0.495, 
0.82, and 0.745, respectively. High Cooper-pair density is 
given by dark regions and low Cooper-pair density by light 
regions. (d,e) The phase of the order parameter for the mul- 
tivortex states with L = 5 at Ho/Hc2 = 1.27 and with L = 6 
at Ho/Hc2 = 1.345. Phases near zero are given by light re- 
gions, phase near 2-w by dark regions. 



FIG. 8. (ar-d) The Cooper-pair density for the multivortex 

states in a disk with vorticity L = 2, 3, 4, 5 at Hq/Hc2 ~ 0.495, 
0.62, 0.965 and 0.82, respectively. High Cooper-pair density 
is given by dark regions, low by light regions. 



FIG. 9. The magnetic field distribution for the square for 
the state with vorticity L = 2 at Ho/Hc2 = 0.42 (a), 0.52 
(b) and 0.62 (c), and with vorticity L = 3 at Ho/Hc2 = 0.62 
(d), 0.72 (e) and 0.82 (f), respectively. High magnetic field is 
given by dark regions and low by light regions. 



FIG. 10. The magnetic field distribution for the square for 
the state with vorticity L = 4 at Ho/Hc2 = 0.72 (a), 0.82 (b), 
0.92 (c) and 1.02 (d). High (low) magnetic fields are given by 
dark (fight) regions. 

FIG. 11. The magnetic field distribution for the L = 11 

state in a square at Hq/Hc2 ~ 1-9 (a) and 1.95 (b). High 
(low) magnetic fields are given by dark (light) regions. 

FIG. 12. The magnetic field distribution for the L = 4 
state in a square for diff'crcnt values of z; z/^ = 0.0 (a), 0.1 
(b), 0.3 (c), 0.6 (d), 1.0 (e) and 10.0 (f), respectively. The 
applied magnetic field is Ho/Hc2 = 0.77. High magnetic field 
is given by dark regions, low magnetic field by light regions. 



FIG. 13. (a-d) Vector plots of the supercurrent in the su- 
perconducting square and (e-h) the contour plots of the phase 

of the order parameter for the L = 1 state at H{)/Hc2 = 0.27 
(a,e), the L = 2 state at Ho/Hc2 = 0.42 (b,f), the L = 3 
state at Ho/Hc2 ~ 0.67 (c,g) and the L = 4 state at 
Ho/Hc2 = 0.745 (d,h). Phases near 2-k are given by dark 
regions, phases near zero by light regions. 

FIG. 14. (a-c) Vector plots of the supercurrent in the super- 
conducting triangle and (d-f) the contour plots of the phase 
of the order parameter for the L = 1 state at Ho/ Hc2 = 0.27 
(a,d), the L = 2 state at Ho/Hc2 = 0.495 (b,e), and the L = 3 
state at Ho/Hc2 = 0.82 (c,f). Phases near 2-it are given by 
dark regions, phases near zero by light regions. 

FIG. 15. The magnetic field range over which the vortex 
states with vorticity L are stable as a function of the vortic- 
ity L for the disk (open circles, solid curves), the square (open 
squares, dashed curves) and the triangle (open triangles, dot- 
ted curves). 

FIG. 16. (a) The expulsion field and (b) the penetration 
field as a function of the vorticity L for disk (solid curves, 
circles), the square (dashed curves, squares) and the triangle 
(dotted curves, triangles). 

FIG. 17. The H - T phase diagram for the disk 
(solid curve), the square (dashed curve) and the triangle 
(dash-dotted curve). Only the superconducting/normal tran- 
sition Hc3 is shown as a function of temperature. The black 
dots indicate the transitions. The inset shows the H — T phase 
diagram for the states with vorticity L = and L — 1. The 
thick curves are the superconducting/normal transitions and 
the thinner curves indicate the expulsion and the penetration. 
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